Entanglement in spatially inhomogeneous many-fermion systems 
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We investigate entanglement of strongly interacting fermions in spatially inhomogeneous environ- 
ments. To quantify entanglement in the presence of spatial inhomogeneity, we propose a local-density 
approximation (LDA) to the entanglement entropy, and a nested LDA scheme to evaluate the entan- 
glement entropy on inhomogeneous density profiles. These ideas are applied to models of electrons 
in superlattice structures with different modulation patterns, electrons in a metallic wire in the pres- 
ence of impurities, and phase-separated states in harmonically confined many-fermion systems, such 
as electrons in quantum dots and atoms in optical traps. We find that the entanglement entropy of 
inhomogeneous systems is strikingly different from that of homogeneous systems. 
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An intense interdisciplinary research effort is currently 
directed at analyzing entanglement in a wide variety of 
physical systems. Apart from the intrinsic interest in 
entanglement as one of the most counterintuitive predic- 
tions of quantum mechanics, this research effort is largely 
motivated by the prospect of simulating or exploiting en- 
tanglement as a resource for quantum information pro- 
cessing and computing [![. 

Particularly promising systems for quantum compu- 
tation are solid-state devices, because of their scalabil- 
ity and the possible integration with existing silicon- 
based technology. However, real solids and solid-state 
devices are necessarily inhomogeneous, i.e., characterized 
by boundaries, interfaces, spatial modulations of system 
parameters, impurities, defects, and externally applied 
fields (arising, e.g., from gate electrodes or from spatially 
confining potentials), among others. Systems of optically 
trapped atoms are also necessarily inhomogeneous, due 
to the trapping potential. If quantum information pro- 
cessing or computing is ever to become reality outside the 
laboratory, we must be able to quantify entanglement in 
realistic, spatially inhomogeneous situations. Unfortu- 
nately, a simple and reliable prescription to quantify and 
interpret the degree of entanglement in inhomogeneous 
systems is still missing. Here we propose a solution to 
this problem in the context of the entanglement entropy. 

The entanglement entropy is defined for a sys- 
tem divided in two subsystems, A and B, as £ = 
Tr a (pa log 2 Pa), where pa = Trgp is the reduced den- 
sity matrix of subsystem A, and p is the density matrix of 
the full system. The Hohenberg-Kohn theorem [2|, y, |4| 
guarantees that £ is a functional of the ground-state den- 
sity, £[n(af)], where x represents sites or spatial coordi- 
nates and n(x) the spatial distribution of particles over 
these sites, but that functional is not known in general. 
In the following, we make the realistic assumption that 
p or £ can be calculated, at least approximately, for a 
spatially uniform system, such as a lattice in which all 
sites are equivalent, or a uniform continuum system with- 



out boundaries. In such spatially homogeneous systems, 
the density n reduces to a constant, and the functional 
£[n(2;)] becomes the function £ hom (n). Our aim is to 
develop a simple, yet reliable, approximation scheme for 
obtaining the ground-state entanglement entropy of in- 
homogeneous systems, £[n(:r)]. 

To this end we propose the local-density approximation 

£[n(x)] « £ LDA [n{x)] = J dx£ hom {n) n _ n(x)l (1) 

in which the functional £[n(x)] is approximated by eval- 
uating the per-site entropy of the homogeneous system, 
£hom ^ a £ foe density distribution of the inhomogeneous 
system. This local-density approximation (LDA) for 
the entanglement entropy is inspired by the LDA made 
in density-functional theory (DFT) for the exchange- 
correlation energy 0, 0| • There is also a conceptual rela- 
tion to the LDA for the thermodynamic entropy made in 
classical DFT of inhomogeneous liquids || and in DFT 
for thermal ensembles @ . 

Clearly, the LDA for the entanglement entropy is a 
general concept, applicable to a wide variety of systems, 
but to be specific, and to compare to previous work, we 
here consider interacting fermions on a lattice, described 
by the Hubbard model, 

H = -tyXcl^Cj+x^ +H.c.) +U ^h^n tl + ^v lr7 n irT 

i,a i icr 

(2) 

where hi a = c\ a Ci a , with a =\, j, is the local density op- 
erator expressed in terms of fermionic creation and anni- 
hilation operators, U is the on-site interaction, t the inter- 
site hopping (below taken to be the unit of energy), and 
Vi a represents external electric and magnetic fields. In- 
clusion of this last term makes the model inhomogeneous. 
The HK theorem guarantees that the global density dis- 
tribution n = rii, ■■! n L determines the ground-state wave 
function of the inhomogeneous L-site model. 

The homogeneous Hubbard model (vi C = 0) accounts, 
approximately, for spin and charge degrees of freedom, 
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FIG. 1: Panel (a): Site-resolved entanglement entropy £ hom (rii, U) = £ m (n, t/)| n _» ni of an open Hubbard chain with impurity 
potential vi — — 1 on the central site. Crosses: exact data for £(rii,U), obtained by fully numerical diagonalization of the 
Hamiltonian. Circles: £ horn (rii,U) evaluated on exact densities. Squares: £ hom (ni,U) evaluated on BA-LDA densities. Panel 
(b): Average per-site entropy £ LDA [rii, U]/L of two superlattices, evaluated on BA-LDA densities. Circles: constant repulsive 
on-site potential v = 3 and repulsive on-site interaction modulated with amplitude A. Squares: constant repulsive on-site 
interaction U — 3 and on-site potential modulated with amplitude A. In both cases the system has L — 200 sites, N = 100 
particles, and the periodicity of the modulation is AL — 20 lattice sites. Panel (c): Average per-site entropy £ LDA [rii, U]/L 
of itinerant interacting particles in the presence of a localized impurity of strength vj. All calculations were done for open 
boundary conditions (OBC), but results for periodic boundary conditions are qualitatively the same. 



for the itineracy of the charge carriers, and for their inter- 
actions. Its rich phase digram and nontrivial physics, to- 
gether with the available exact solution in one dimension, 
has attracted much interest also in the quantum infor- 
mation community. In particular, several groups inves- 
tigated entanglement measures for the Hubbard model 
(see, e.g., Refs. @, S, @, EU EL E El and references 
therein). As a result, the single-site entanglement en- 
tropy (in which subsystem A is one site) of the homoge- 
neous Hubbard model in one dimension is known analyt- 
ically 0, S, LH : 



£ hom (n, U) -- 

2+ m log2 
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.(3) 



Here £ hom (n,U) is the single-site entanglement entropy 
(Shannon- von Neumann entropy of four probabilities) 
per site, and e = E(n,U)/L is the ground-state energy 
per site, both as functions of interaction strength U and 
particle density n = N/L, where N and L are the total 
number of particles and lattice sites, respectively. Since 
e(n, U) can be obtained [3, EH from the Bethe-Ansatz 
integral equations, Eq. (J3J) can be used to extract detailed 
quantitative information about entanglement in various 



phases and at the transitions between them 0, H, [Hj]. 
However, this expression is valid only for the homoge- 
neous Hubbard model, i.e., one in which all sites are 
equivalent. Such idealization is most useful for analytical 
and numerical work, but provides only a rather imperfect 
picture of the physical situation in real solids or devices. 

To make use of expression ([3]) for the entanglement en- 
tropy of homogeneous Hubbard models also in inhomo- 
geneous situations, we apply the local-density approxi- 
mation {1} in the form 



£[n h U] w £ LDA [n t , U] = £ £ hom (n, U)\ r 



(4) 



For a given distribution m of particles over sites, 
Eq. Q can be evaluated immediately. If, on the other 
hand, this distribution is not known, it must be obtained 
in a separate calculation. This is the case, e.g., if the 
system is specified by giving the potential the particles 
move in, instead of the spatial distribution of particles. 
Experimentally, specifying the potential is more realistic, 
and we thus focus on this case. 

To obtain the ground-state charge and spin densities 
for a given external potential, we again appeal to density- 
functional theory, this time employing the LDA in its 
standard formulation 0,0], i.e., as an approximation to 
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the exchange-correlation energy that is to be added to 
the mean-field energy functional. Minimization with re- 
spect to n(x) then yields the ground-state energy and 
ground-state density profile. Specifically for the Hub- 
bard model, we use the Bethe-Ansatz LDA (BA-LDA) of 



Refs. [1J,|15|]. The combined calculation is thus a nested 



LDA, where BA-LDA is used to generate selfconsistent 
density profiles, and the entropy LDA fl} is used to pre- 
dict the entropy on the resulting inhomogeneous density 
distribution. 

Both BA-LDA and the entropy LDA by construction 
become exact for a spatially uniform system. BA-LDA 
has been shown by comparison to Monte Carlo, Bethe 
Ansatz, and density-matrix renormalization group data 
to provide density profiles that agree within a few per- 
cent with those of other many-body methods, even for 
systems far from the uniform limit 14 , T^, • To quan- 
tify the reliability of the entropy LDA, we turn to small 
finite chains, which are far from the uniform limit, but 
for which the exact density matrix p can be obtained 
by full numerical diagonalization of the Hamiltonian, so 
that the exact entanglement entropy £ is known. Since 
small systems are a worst-case scenario for LDAs (which 
derive from the thermodynamic limit), this is a partic- 
ularly severe test for the entropy LDA and nested LDA 
concepts. 

Panel (a) of Fig. [T] displays the site-resolved entan- 
glement entropy of an open chain chain with an impu- 
rity at the central site. Exact data are obtained by nu- 
merically diagonalizing the many-body Hamiltonian @ 
and constructing the density matrix and the entangle- 
ment entropy from the exact eigenfunctions. This pro- 
cess also yields the exact density profile, which we use to 
evaluate the function £ hom (n, U) entering the entropy 
LDA, Eq. 0$. The deviation between both sets of data 
measures the quality of the entropy LDA. In a separate 
calculation, we also generate BA-LDA density profiles 
3 [H EH, and evaluate £ hom (n, U) on them. This re- 
sults in the nested LDA procedure, which can also be 
applied where the exact solution is not available. As 
the figure shows, all three approaches agree closely, too 
within 1%, even far away from the limit in which LDAs 
become exact. 

Having established the viability and reliability of the 
LDA for the entanglement entropy and of the nested 
LDA, we now turn to further experimentally relevant in- 
homogeneous systems. First, we investigate the effect a 
spatial modulation of system properties, in the form of a 
superlattice, has on the entanglement. Naturally occur- 
ing or man-made systems displaying spatial modulations 
of their properties on a nanoscale are important both as 
paradigms of simple nanotcchnological devices, and as 
particular examples of emerging spatial inhomogeneity 
in strongly interacting many-electron systems. To model 
superlattice structures in the Hubbard model we follow 
Ref. [161 ], which, however, did consider only energies, not 



entropies. 

Panel (b) of Fig. [1] displays the behaviour of £ LDA [rii] 
for two representative superlattice structures. According 
to the nested LDA concept, the data were obtained by 
first running an ordinary BA-LDA calculation for a fixed 
distribution of on-site interaction and potentials, in order 
to generate the densities ni. In a second step the resulting 
self-consistent density profile is substituted in Eq. ([I]) , in 
order to obtain £ LDA [ni], for the specified values of Ui 
and Vi. 

Two things leap to the eye in Fig.QJb): First, the larger 
the amplitude of the modulation, the less entangled the 
system becomes. Indeed, the spatial inhomogeneity of 
the superlattice is expected to disrupt the entanglement 
present between itinerant charge carriers in the homoge- 
neous system. Second, modulations in the on-site poten- 
tial have a much more drastic effect on the entanglement 
entropy than those of the on-site interaction. Clearly, lo- 
cal fluctuations in the potential disentangle the system 
much more efficiently than local fluctuations in the inter- 
action. This is consistent with what was previously [l6| 
observed for the ground-state energy, and again points to 
the importance of local electric fields in strongly corre- 
lated systems [l7[ . It also means that local electric fields 
are rather effective in reducing the degree of entangle- 
ment between itinerant charge carriers. 

Next, we turn to impurity systems. Specifically, we 
inquire what effect a local impurity potential has on the 
entanglement between itinerant electrons. We model the 
impurity as Vi = vjSi^/^}, i-e., place a local electric po- 
tential on the central site of the L-site chain. Panel (c) 
of Fig. [T] shows that both for attractive (vj < 0) and 
repulsive (vi > 0) impurities, entanglement between the 
itinerant particles is diminished, but by a much smaller 
margin than in the superlattice case, where the potential 
is periodically repeated instead of acting on just one site. 
Stronger on-site interactions U tend to freeze the density 
distribution and thus reduce the number of degrees of 
freedom, which lowers the entanglement entropy, as dis- 
cussed for homogeneous systems in Ref. 13j. For small 
U, attractive impurities are more efficient in disrupting 
itinerant entanglement than repulsive ones, but this dif- 
ference disappears for larger U. For attractive impurities, 
we observe a flat structure at small negative values of Vi . 
For these values, the average density of the impurity site 
remains fixed at m = 1, which means that the strenght of 
the impurity potential (which couples to the density) be- 
comes irrelevant, until it increases enough to draw more 
than one particle to the impurity site, at which point the 
frozen site again becomes available as a degree of freedom 
and £ is slightly increased. This freezing of the density 
distribution at n = 1 is a strong-interaction effect (it van- 
ishes for U = 0), similar to the Mott insulator in solids, 
and to Coulomb blockade in quantum dots [3, , and 
to the density plateaus in our next example. 

Finally, we consider harmonically confined fermions, 
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FIG. 2: Thick curve: Average per-site entropy £ LDA (rn, U)/L 
evaluated on BA-LDA densities, of harmonically trapped 
fermions on a lattice, as a function of curvature of the confin- 
ing potential. Thin curve: numerical derivative of the thick 
curve. Insets: density profiles corresponding to special points 
on the two curves. 

such as electrons in a quantum dot [l!^ or fermionic 
atoms in an optical trap [2(| • Figure [2] displays the en- 
tanglement entropy and its derivative as a function of 
curvature of the trapping potential. Increasing confine- 
ment diminishes entanglement of the confined particles, 
but not in a uniform way. As the curvature becomes 
larger, the system passes through several states with dis- 
tinct density profilesjsee insets), whose nature has been 
clarified elsewhere (la . Ell ]. Interestingly, each transi- 
tion corresponds to distinctive features in the entangle- 
ment entropy, many of which are more clearly visible in 
its derivative: the switching on of the confining poten- 
tial leads to an initial drop in entanglement, a spike in 
the derivative signals the appearance of a Mott insulat- 
ing region in the trap center, a smaller peak indicates 
ressurection of a metal-like region in the center, and an- 
other broad peak is associated with formation of a band- 
insulating state at the trap center. The entanglement 
entropy can thus be used as a tool for characterization 
of the various possible states, even though these are not 
phases in the thermodynamic sense. 

On the other hand, it also becomes clear that the ac- 
tual value of E , and its variation with system parame- 
ters, depend in a highly nontrivial way on fine details of 
the system. This dependence is a severe obstacle for at- 
tempts to quantify and exploit entanglement in real-life 
environments. Note that this issue is conceptually dis- 
tinct from the more commonly discussed problems due 



to decoherence, which also occurs in spatially homoge- 
neous situations. 

In summary, we find that for all types of inhomogeneity 
investigated here — superlattice structures with differ- 
ent modulation patterns, itinerant electrons in a metallic 
wire in the presence of impurities, and phase-separated 
states in harmonically confined many-fermion systems — 
the entanglement entropy of inhomogeneous systems is 
strikingly different from that of homogeneous systems. 
The influence of spatial inhomogeneity is an essential as- 
pect of entanglement in real systems, which cannot be 
simulated or understood by approaches based exclusively 
on uniform systems. This influence must be understood 
and quantified if entanglement is to be used as a resource 
for performing quantum information processing in actual 
devices. The entropy LDA and the nested LDA concept, 
proposed here, may be useful tools for this investigation. 
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